Analysing solar-like oscillations with an automatic pipeline 
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Abstract. Kepler mission will provide a huge amount of asteroseismic data during the next few years, among which 
hundreds of solar-like stars will be targeted. The amount of stars and their observation length represent a step forward in 
the comprehension of the stellar evolution that has already been initiated by CoRoT and MOST missions. Up to now, the slow 
cadence of observed targets allowed an individual and personalized analysis of each star. During the survey phase of Kepler, 
this will be impossible. This is the reason why, within the AsteroFLAG team, we have been developing automatic pipelines 
for the Kepler solar-like oscillation stars. Our code starts by finding the frequency-range where p-mode power is present and, 
after fitting the background, it looks for the mode amplitudes as well as the central frequency of the p-mode hump. A good 
estimation of the large separation can thus be inferred in this region. If the signal to noise is high enough, the code obtains the 
characteristics of the p modes by doing a global fitting on the power spectrum. Here, we will first describe a few features of 
this pipeline and its application to AsteroFLAG synthetic data to check the validity of the code. 
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INTRODUCTION 

Kepler mission, which is dedicated to the search for exo- 
planets as well to asteroseismic analysis was launched 
in March 2009. During the survey phase, a few hun- 
dreds of solar-like oscillation stars will be observed 
for one month. To analyse this huge number of stars 
very quickly, several pipelines have been developed. We 
present here our pipeline applied to stars simulated by 
the AsteroFLAG group. It has also been tested on a few 
solar-like oscillation stars already observed with CoRoT, 
such as HD49933 0, HD181906 0, and HD181420 
0. 

PIPELINE DESCRIPTION 

The pipeline is constituted of five main packages: 

1. P-Mode Range Search: we look for the frequency 
range of the p-mode excess power as well as the 
large-frequency separation estimation: 

• first, we estimate the large separation, Av, by 
calculating the Power Spectrum (PS) of PS 
from 200 to 6000 jliHz (called PS2 hereafter, 
see an example on Fig. 3). 

• then, we take sub-regions of 600 /xHz in the 
PS, compute the PS2 and look for the relative 



power of the highest peak close to the value of 
Av found. By calculating the threshold corre- 
sponding to the probability of having a peak 
at this position, we obtain the range of p-mode 
excess power (/„„■„ and f max ) with a given con- 
fidence level (varying from 70 to 95% 0) 

2. Background Fitting: we fit the background with 3 
components, 2 harvey laws and a flat photon 
noise. The technique used is the maximum likeli- 
hood minimization. 

3. Bolometric amplitude of the £ = mode: we look 
for the amplitude envelope of the PS 

• we correct the PSD from background and we 
smooth it over a window of 1 or 2xAv; 

■ we fit the result with a Gaussian to get P mux and 

• we convert the power to the amplitude per ra- 
dial mode A max (6) and to bolometric ampli- 
tude 0. 

4. First determination of the Mass and Radius of the 
star: we use the inversed scaling laws [8, 9| 

5. Global fitting: we determine the table of guesses 
of modes frequencies from v majc and Av. We use a 
maximum likelihood minimization fitting of all the 
modes in a given frequency range including also the 
background. The code allows several possibilities 



like including Bayesian constraints, different mode 
identification and, I = 3 modes. 



The pipeline also contains three more packages: 

1. Curvelets: Enhance the ridges of the p-mode signal 
in the echelle diagram ifTUl 

2. PSSPS (Power Spectrum of a Short part of the 
Power Spectrum), which is also part of package #1: 

• we estimate Av from the standard stellar pa- 
rameters 

• we search for the spacing in the PS of a short 
slice of the PS where the p-mode excess power 
appears ifTTI 

3. Determination of the rotation period (with several 
techniques among which the wavelets) 

RESULTS OBTAINED WITH THE 
ASTEROFLAG HARE AND HOUNDS 

For 3 different stars generated by the asteroFLAG team, 
we have applied our pipeline (from packages #1 to 4): 
Pancho, Boris and Arthur [ 12 1. Each star has been simu- 
lated with 3 different magnitudes (m,,) (or Signal to Noise 
Ratio, SNR): 9, 10 and 1 1. They are all dwarves (V). Pan- 
cho has a r e /y=6383 K, for Boris, r e yy=5780 K and for 
Arthur, 7^=6420 K. Each time-series is 3 years long. In 
our analysis, we used sub-series of 3 1 days as KEPLER 
will observe this length of data during the survey phase 
before focusing on a few star for several years. 

Fig. 1 shows the results for the estimation of Av with 
package # 1 for magnitudes of 9 and 11. We can see that 
in most of the cases, we manage to find a good estimation 
of Av. For m v =9, in 97% of the cases, we obtained a 
good value with uncertainties lower than 3 flHz. For 
m v =10, a good estimation is obtained for 65% of the 
cases and 36% for m,.= ll. Once we have an estimation 
on the large separation, we calculate the frequency-range 
of the p-mode excess power. Table 1 gives a statistical 
overview of the success ratio concerning the detection of 
the frequency range where most of the p-mode bump is 
found. For a bright star, almost 100% of the power excess 
is correctly found. But above a SNR of 11, results start 
to become less and less reliable. 

TABLE 1. Percentage of good results 
for the frequency-range of the p-mode 
excess power 
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FIGURE 1. Histograms of the results of the pipeline on the 
estimation of Av. For each m v , we have used 3 stars. 



Finally, after removing the background fit for each 
star, we calculate the maximum Amplitude of 1=0 mode. 
Here, we took only the stars where the frequency range 
was correctly found. The results for m v =9 and 11 are 
shown on the histograms of Fig. 2. The amplitude is also 
well retrieved for most of the stars with SNR of 9 and 10. 
The error bars associated vary from 0.8 to 2ppm for the 3 
SNR. Fig. 3 and 4 show the results obtained for 31 days 
of Pancho and m v = 9. 

CONCLUSIONS 

The pipeline applied to the asteroFLAG hare and hounds 
datasets gives good results concerning stars with SN of 9 
and 10. Within the error bars we retrieve the values of the 
main global parameters of the modes. This pipeline for 
the packages 1 to 4 is now working in an automated way 
and was tested on CoRoT targets as well as on simulated 
stars from the Kepler-asteroFLAG simulator. 
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FIGURE 2. Histograms of the results of the pipeline on the 
maximum amplitude per radial mode. Same as figure 1. 
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FIGURE 3. Example of the PS2 for 3 1 days of Pancho. 
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FIGURE 4. Example of Pancho: results of package #1 (top) 
and 4 (bottom). 
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